Astronomy & Astrophysics manuscript no. ttau_v39 ©ESO 2015 

May 26. 2015 


Letter to the Editor 


Examining the T Tauri system with SPHERE * 

Gergely Csepany 1,2 , Mario van den Ancker 1 , Peter Abraham 2 , Wolfgang Brandner 3 , and Felix Hormuth 3 

1 European Southern Observatory, Karl-Schwarzschild-Str. 2, 85748 Garching bei Miinchen, Germany 

2 Konkoly Observatory of the Research Centre for Astronomy and Earth Sciences, Hungarian Academy of Sciences, Konkoly Thege 
Miklos ut 15-17, 1121 Budapest, Hungary 

3 Max-Planck-Institut fiir Astronomie, Konigstuhl 17, 69117 Heidelberg, Germany 

Received / Accepted 


in 

t-H 

o 

(N 

03 


in 

(N 

P? 

cn 

Ph 

6 

C/3 


ABSTRACT 

Context. The prototypical low-mass young stellar object, T Tauri, is a well-studied multiple system with at least three components. 
Aims. We aim to explore the T Tau system with the highest spatial resolution, study the time evolution of the known components, and 
re-determine the orbital parameters of the stars. 

Methods. Near-infrared classical imaging and integral field spectrograph observations were obtained during the Science Verification 
of SPHERE, the new high-contrast imaging facility at the VLT. The obtained FWHM of the primary star varies between 0.050" and 
0.059", making these the highest spatial resolution near-infrared images of the T Tauri system obtained to date. 

Results. Our near-infrared images confirm the presence of extended emission south of T Tau Sa, reported in the literature. New 
narrow-band images show, for the first time, that this feature shows strong emission in both the Br-y and H 2 1-0 S(l) lines. Broadband 
imaging at 2.27 pm shows that T Tau Sa is 0.92 mag brighter than T Tau Sb, which is in contrast to observations from Jan. 2014 
(when T Tau Sa was fainter than Sb), and demonstrates that T Tau Sa has entered a new period of high variability. The newly obtained 
astrometric positions of T Tau Sa and Sb agree with orbital fits from previous works. The orbit of T Tau S (the center of gravity of Sa 
and Sb) around T Tau N is poorly constrained by the available observations and can be fit with a range of orbits ranging from a nearly 
circular orbit with a period of 475 years to highly eccentric orbits with periods up to 2.7 x 10 4 years. We also detected a feature south 
of T Tau N, at a distance of 144 ± 3 mas, which shows the properties of a new companion. 
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T Tauri (HD 284419), initially considered to be the prototype of 
a class of low-mass pre-main sequence stars, is now known to 
be a complex young system of at least three stellar components 
(T Tau N, and the close binary system T Tau Sa/Sb) as well as of 
jets, and outflows. Whereas T Tau N suffers little extinction ( A v 
= 1 m 46; White & Ghez|200l l, T Tau Sb is hidden behind about 
15 mag of extinction, while T Tau Sa appears to be even more 
obscured (van den Ancker et al.||1999] jvan Boekel et al.||2010j 


presence of a circumstellar disk viewed nearly edge-on around 
this source as well (Duchene et al. 2005). In the high-resolution 


images van Boekel et al.|(|2010|) found circumstellar emission 


around T Tau Sa in the K -band on scales of 5 to 20-30 AU. It is 
as of yet unclear whether this emission can be directly associated 
with the tentative edge-on disk of T Tau Sa. All three compo¬ 
nents of the T Tau system are actively accreting and are spatially 


|Duchene et al.|2005|. 

Orbital motion of Sa and Sb has been measured over the 
last 16 years, yielding a well-constrained orbital model for Sa- 
Sb with a semi-major axis of 89 mas and a period of 29 years 
( [Schaefer et al. 2014| >. Using the orbital data, the mass of 
T Tau Sa was shown to be 1.9-2.3 M 0 and the mass of Sb is be¬ 
tween 0.7 and 0.9 M Q ( |Duchene et al.|2006l|Kohler et al.|2 008 ). 
The orbital motion of the Sa-Sb system around T Tau N is less 
well constrained, with initial fits yielding an orbital periods be¬ 
tween 475 and 23000 years, and a total system mass between 4.5 
and 5.6 M Q ^Johnston et al.|2004[|Kohler et al.|2008| l. 

Millimeter interferometry first revealed the presence of a 
circumstellar disk around the optically dominant component 
T Tau N ( |Hogerheijde et al.|1997[[Akeson et al. 1998| l. The pres¬ 
ence of CO absorption in the spectrum of T Tau Sa suggests the 


unresolved hydrogen recombination line emitters (Kasper et al. 
|2002[|Gustafsson et al.|20K)] i. 

In this paper we present new infrared images and spec¬ 
troscopy of the T Tau system, obtained with SPHERE at the 
VLT. At small scales, these new images present the highest con¬ 
trast data obtained so far of T Tau. We discuss the orbital and 
photometric evolution of the known components of T Tau and 
present evidence for the presence of a new tentative companion 
south of T Tau N. 


2. Observations and Data Analysis 

Observations of the T Tau system were acquired during the 
Science Verification of SPHERE0 the newly installed extreme 
Adaptive Optics facility at the VLT. The first set of observations 
was obtained on Dec. 9 th , 2014, when SPHERE was operated in 
the IRDIS classical imaging mode with the ContJ (1.216 pm). 
Pa -p (1.281 pm), H 2 (2.128 pm), Br-y (2.163 Atm) and CntK2 


* Based on observations made at the La Silla Paranal Observatory 
under programme ID 60.A-9363(A) and 60.A-9364(A). 


http://www.eso.org/sci/facilities/paranal/instruments/sphere.html 
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Fig. 1. Classical images of the T Tauri system in different filters 
with logarithmic intensity scaling. From left top to bottom right: Pa-/3 
(1.281 pm), Hi (2.128 pm), Br-y (2.163 pm), and CntK2 (2.270 pm). 


(2.270 pm ) filters. Calibration frames were acquired for the ge¬ 
ometrical distortion map, which were used during the data pro¬ 
cessing to ensure astrometric precision. A second set of obser¬ 
vations was acquired on Jan. 22 nd , 2015, when SPHERE was 
operated in the IRDIFS mode with a 155 milli-arcsecond (mas) 
diameter apodized Lyot coronagraph, offering simultaneous In¬ 
tegral Field Spectrograph (IFS) observations from 0.95-1.65 pm 
and imaging in the Ks (2.181 pm) filter. Observations with the 
star offset from the coronagraphic mask were also obtained in 
order to be able to measure total fluxes. A PSF reference star 
(TYC 1290-457-1) was observed right after the IRDIFS obser¬ 
vation of T Tau, at similar airmass and atmospheric conditions 
as the science target. 

The obtained data sets were reduced with the pre-release ver¬ 
sion 0.14.0-2 of the SPHERE pipeline. Apart from the usual 
steps of bias subtraction, flat-fielding and wavelength calibra¬ 
tion, the pipeline also includes a correction for geometric dis¬ 
tortions. Post-processing of the data was done using the FITSH 
software package ( |Pal|2012 i to sum the data resulting from the 
two IRDIS channels. The achieved FWHM of the diffraction- 
limited core of T Tau N varies between 0.050" and 0.059" (from 
ContJ to CntK2), making these the highest spatial resolution 
near-infrared images of the T Tau system obtained to date. The 
classical images are shown in Fig. |T| The H 2 , Br-y and CntK2 
filters suffer from a number of known ghosts which can be seen 
close to bright stars. Figure [2]features the IFS images taken with 
the coronagraph, in which a cross-like pattern can be seen that is 
caused by reflections from the wires on which the coronagraphic 
mask is suspended. The longer wavelength frames in the IFS ob¬ 
servation suffer from a flat fielding error pattern; either the flat 
fields were inadequate or the pre-release version of the pipeline 
failed to properly reduce the data. We omit performing precise 
astrometry and photometry in those frames, and only use them 
to determine the spectrum with a larger aperture to compensate 
for the error pattern. 

Positions and fluxes (both relative to T Tau N) of the stellar 
components T Tau Sa and Sb were measured and are listed in Ta- 
ble[l] Relative fluxes in this Table were computed using aperture 
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Table 1. Positions and magnitude differences (relative to T Tauri N) of 
the stellar components. 


Object 

Filter 

Sep. ( 

mas) 

PA (°) 


Magn. diff. 

T Tau Sa 

h 2 

688 

± 

3 

193.4 

± 

0.2 

2.41 

± 

0.01 


Br-y 

686 

± 

3 

193.4 

± 

0.2 

2.26 

± 

0.01 


CntK2 

685 

+ 

3 

193.3 

+ 

0.2 

1.95 

± 

0.01 


Ks 

687 

± 

15 

193.1 

± 

1.0 

1.82 

± 

0.03 

T Tau Sb 

ContJ 

589 

± 

6 

198.2 

± 

0.4 

6.14 

± 

0.12 


Pa-/? 

590 

± 

6 

198.2 

± 

0.4 

5.77 

± 

0.09 


h 2 

595 

± 

3 

197.9 

± 

0.3 

3.15 

± 

0.01 


Br-y 

593 

+ 

3 

198.0 

+ 

0.3 

3.09 

± 

0.02 


CntK2 

592 

± 

3 

197.9 

± 

0.3 

2.87 

± 

0.02 


Ks 

589 

+ 

15 

197.7 

+ 

1.3 

2.58 

± 

0.05 

Tentative 

1.02 pm 

140.9 

± 

2.7 

199.0 

± 

2.8 

4.96 

± 

0.30 

companion 

1.19 pm 

143.7 

± 

2.5 

198.0 

± 

1.6 

4.39 

± 

0.19 

(IFS obs.) 

1.34 pm 

145.8 

± 

3.1 

197.2 

± 

2.0 

4.37 

± 

0.30 

photometry, with the aperture radius taken equal to the radius 


of the Airy disc (73.5 mas). The listed uncertainties of the stel¬ 
lar positions and angles include the error of the centroid fitting, 
the pixel-scale conversion uncertainty and relative atmospheric 
diffraction. The latter error contribution is calculated by assum¬ 
ing the maximum relative shift in position of the components 
diffraction in the used filter. 


3. Results 


3. 1. Classical Imaging 

The two shortest wavelength images (ContJ and Pa-/?) show 
the primary (T Tau N) and one of the southern components 
(T Tau Sb). In addition, some extended emission may be present 
at low intensities. In the longer wavelength images (H 2 , Br-y 
and CntK2), the second known southern component (T Tau Sa) 
is also visible. In all images, the stellar sources T Tau N, Sa and 
Sb appear indistinguishable from a point-source. 

Observational 3<x detection limits in the different filters were 
calculated from our data and are shown in Fig. [4] (online mate¬ 
rial) by calculating the magnitude of a hypothetical source that is 
hidden in the background noise at the given location. The bumps 
in the plots at ~ 250 and ~ 600 mas are the result of the Airy 
rings. We do not detect the tentative northern component, de¬ 
noted as T Tau (O’) by Nisenson et al. (1985) (and later also 
observed by Maihara & Kataza 1991) . Our upper limits show 
that this tentative component, if exists, must be at least 4.7 mag¬ 
nitudes fainter than T Tau N over the entire 1.2-2.2 pm range 
covered by our observations. 

In our CntK2 (2.27 pm) image, T Tau Sa appears distinctly 
brighter (0.92 mag. Table [TJ than T Tau Sb. This is in marked 
contrast to the latest AT-band brightness ratio measurements, ob¬ 
tained in January 2014, in which Sa was about 0.8 magnitudes 
fainter than Sb ( | Schaefer et al.|20f4) . This variability of T Tau Sa 
of 1.7 magnitudes in less than a year is similar to the photometric 
behaviour shown prior to 2006 ( Beck et al.|2004| Duchene et al. 
2005), and demonstrates that the period of lesser variability re- 


i Boekel et al. {2010 , Schaefer et al. 


ported from 2006-2013 
12014) 1 has come to an enc. 

In addition to the continuum filters, our data set includes 
measurements in several narrow-band filters centered on well- 
known emission lines (Pa-/?, Br-y and the H 2 1-0 S(l) line). 
From the brightness ratio of Br-y to the nearby CntK2 filter we 
conclude that T Tau N (brightness ratio 1.52) shows strong Br-y 
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emission, most likely linked to the accretion onto the central star. 
Measurements for Sa and Sb (brightness ratios of 1.14 and 1.22 
respectively) are inconclusive, although we note that Br-y emis¬ 
sion has previously been detected from both these components 
( jGustafsson et al.|2010| l. 

For Ht, the brightness ratio H2/CntK2 is very high (1.90) for 
T Tau N, demonstrating the presence of strong Hj 1-0 S(l) emis¬ 
sion close to this component. Previous studies found H 2 emis¬ 
sion around T Tau N and S ( |Kasper et al.|2002||Beck et al.|2004[ 
Gustafsson et al. 2010), which our observations also confirm. 
Our CntK2, Br-y and H 2 images (Fig. [5] in the online material) 
show extended emission at a peak flux of 6<x 
level south of T Tau Sa, first described by 
( |2010| Fig. [5j. The ratio of the flux in the H 2 
that in CntK2 (1.78 and 1.53, respectively) is high, demonstrat¬ 
ing that this feature possesses emission in the Br-y and the H 2 
1-0 S(l) lines. 


above background 
van Boekel et al. 
and Br-y filters to 


3.2. Integral Field Spectrograph 


Our IFS data show a point-like feature south of T Tau N, which 
we identify as a tentative new component in the system. This 
new companion is shown in Fig. [2] which shows the PSF ref¬ 
erence star TYC 1290-457-1, the T Tau system and the T Tau 
system with the PSF reference star subtracted at two different 
wavelengths (1.017 and 1.341 pm). The astrometric and photo¬ 
metric measurements of the new companion are computed with 
a 0.1 pm binning, and are listed in Table [I] 

Although several artefacts can be seen in the IFS data shown 
in Fig. [2] the feature we label here as a tentative new compan¬ 
ion is unique in several aspects: (1) the tentative companion is 
only seen in the T Tau IFS data; no similar feature is present in 
the PSF data, (2) the tentative companion has the same position 
at the different wavelengths in the IFS data cube, excluding the 
possibility that this could be a speckle from the AO system, (3) 
the spectrum of the tentative component is much redder than the 
spectrum of T Tau N (the spectral slope in F v between 0.97-1.30 
pm is 1.3 for the companion and 0.9 for T Tau N), excluding the 
possibility that this could be a reflection from T Tau N in the 
instrument’s optics. 

Close inspection of the IRDIS images in Fig. [T| also reveals 
a slightly extended point-like source just south of T Tau N, al¬ 
though its distance is close to the first Airy ring. As we do 
not have a PSF reference observation for the IRDIS frames, we 
cannot perform precise measurements in those frames and we 
rely on the IFS frames to determine the companion’s properties. 
The tentative companion is not detected in the IRDIS Ks im¬ 
ages and is not found in the Kp- and H-band images made by 
Keck/NIRC2 ( jSchaefer et ak][2f)T4| [Duchene et aL|2006] ). The 
IRDIS images have an observational contrast limit of» 3.3 mag¬ 
nitude at 140 mas, one magnitude less than the brightness of the 
tentative companion, and the Keck/NIRC2 images have a con¬ 
trast limit of ss 1.0 at the same position. Therefore it is likely that 
the tentative companion is lost in the sky background in these 
images. 


The absolute magnitude of the tentative companion can be 
calculated using the 2MASS brightness of T Tau N, and the 
distance modulus derived by Loinard et al. (2005[), yielding 


J maq — 6.0+0.4. Assuming an age of the system of 2 Myr (jWhite 


Baraffe et al. 


& Ghez 2001} , and using the evolutionary tracks of 
( 1998| >, this absolute magnitude is compatible with a stellar or 
sub-stellar companion with mass between 0.06 and 0.15 M Q . 
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Fig. 2. The T Tauri system as seen at two different wavelengths in the 
IFS observation. Top: PSF reference star (TYC 1290-457-1), Middle: 
T Tauri system, Bottom: T Tauri with the PSF reference subtracted. 


3.3. Astrometry 


We compared our new astrometric measurements of the relative 
positions of T Tau N, Sa and Sb with all available astrometry 
of the T Tau system from the literature (Fig. [3}. The literature 
contains (mostly recent) astrometric measurements in which Sa 
and Sb have been spatially resolved, and (mainly older) measure¬ 
ments in which only one southern component was detected. In 
our IFS observations T Tau Sa becomes visible at wavelengths 
longer than 1.4 pm, consistent with the IRDIS images, where 
T Tau Sa is only visible in the longer wavelength observations. 
Therefore we can safely assume that any observation at wave¬ 
lengths shorter than 1.4 pm that did not resolve the southern 
components, has in fact observed T Tau Sb. We thus interpret 
those measurements as reflecting the position of T Tau Sb. We 
also know that the luminosity of T Tau Sa is comparable to T 


Tau Sb in the K-band and highly variable (van Boekel et al. 2010 
Scha efer et al.|2 014). This makes it impossible to use unresolved 
K-band observations to determine the exact position of Sa and 
Sb. We therefore omit from our analysis any literature observa¬ 
tions in the K-band, in which Sa and Sb have not been resolved. 

As can be seen from Fig. [3] significant orbital motion of the 
Sa-Sb binary around T Tau N, and in the Sa-Sb system itself 
has been detected. The positions are relative to either T Tau N 
(left) or T Tau Sa (right). The observations are labelled after the 
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Table 2. New orbital parameters for T Tau N-S. 


Parameter 

Range 

a (") 

0.57 

- 16.8 

a (AU) 

84 

- 2465 

Eccentricity 

0.10 

- 0.98 

Inclination (°) 

0.6 

- 78.1 

wo 

1.2 

- 178.0 

fin 

14.1 

- 176.7 

Period (years) 

475 

- 27000 

Tperiastron (years) 1400 

- 7200 


literature source. The unresolved observations are shown in light 
grey for the K-band, and blue for the H-band. Both are corrected 
for the motion of T Tau Sb around T Tau Sa, using the orbital 
parameters from |Schaefer et al.| ( |2014| ). It can be seen that the 
H-band data align with the resolved observations, whereas the 
K-band data (which are affected by the variable brightness ratio 
of Sa/Sb) seem to deviate. 

The corrected H-band points and the epochs from resolved 
observations miss the orbital fit of the T Tau S around T Tau N 
(the center of gravity of Sa and Sb) from Kohler et al.]( j2008| l, 
clearly showing that this previously derived orbit does not cor¬ 
rectly take into account the orbital motion of T Tau Sb around 
Sa. Therefore we searched for new orbital fits. We have excluded 
some observations from the orbital fit as their uncertainties were 
too low to acquire a xl < 1 fit, as shown in Fig. [ 3 ] There are 
also radio observations of this system, which list astrometric po¬ 
sitions for presumably T Tau Sb, see Johns ton et al.| ( |2003) and 
Smith et al.]( |20031 >. However, even after correcting for the orbital 


motion of Sb around Sa, the radio observations do not align with 
the optical observations, and rather group with the unresolved 
K-band observations. This suggests that the origin of the radio 
emission may be more complex, and it may not originate only 
from T Tau Sb, therefore we omitted them from the plot and the 
orbital fit. 

A random grid search of the parameter space yields a wide 
range of possible solutions for the orbit of T Tau S around T 
Tau N. The possible range of orbital parameters is shown in Ta¬ 
ble 2] We note that some of the parameters listed in this table 
exc ude the orbital solution derived by Kohler et al.|( 2008) ). This 
is a natural consequence of the clear differences in our fits as 
can seen from Figure [3] We note that, assuming the same mea¬ 
surement precision that we achieved with SPHERE, a unique 
solution for the orbit should be possible with new data obtained 
around 2027-28. Regarding the orbital motion of Sa around Sb, 


our fit is essentially identical to the orbital fit by Schaefer et al. 
( f20T4| >. 


4. Summary and Conclusions 

In the new SPHERE data on the T Tauri system presented here, 
two components, T Tau N and Sb are visible in all classical im¬ 
ages, while Sa only becomes visible at wavelengths longer than 
1.4 yum. All components appear indistinguishable from a point- 
source. We detected a new tentative companion south of T Tau 
N, at a position of 144 ± 3 mas, PA: 197 + 2°, relative to T Tau N. 
This companion is not present in the PSF reference observation, 
and its spectrum differ from the other stellar companions in the 
system. If this is a new stellar or sub-stellar source within the 
system, we estimate its mass to be between 0.06 and 0.15 M G . 

New astrometric measurements of the relative positions of 
T Tau N, Sa and Sb were compared with astrometric data from 
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the literature. The astrometric epochs agree with the orbital fit 
from Schaefer et al. ( |2014j ) for the expected orbital motion of 
Sb around Sa. The orbit of the T Tau S system around T Tau N 
is poorly constrained by the existing observations and can be fit 
with a range of orbits ranging from a nearly circular orbit with a 
period of 475 years to highly eccentric orbits with periods up to 
2.7 x 10 4 years. 

T Tau Sa appeared 0.92 magnitude brighter than T Tau Sb in 
the CntK2 image, which is in contrast to the latest K-band bright¬ 
ness ratio measurements from January 2014 (Schaefer et al. 
2014)1, where Sb was brighter than Sa. This photometric variabil 


ity is similar to the behaviour shown prior to 2006 (Beck et al. 
2004; Duchene et al. 2005), and thus demonstrates that the pe¬ 
riod of lesser variability, when Sa and Sb had comparable bright¬ 
nesses between 2006-2013, has come to an end. The brightness 
ratio H2/CntK2 is 1.90, very high for T Tau N, demonstrating 
the presence of strong H 2 1-0 S( 1) emission close to this compo¬ 
nent. Previous studies found H 2 emission around T Tau N and S 
|Gustafsson et al.)(]20iO|, which our observations also confirm. 
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Fig. 3. Left: Epochs of T Tau S around T Tau N. The coordinates are relative to T Tau N. Observations in which the Sa/Sb pair is not resolved 
are shown in either light gray or blue. The coloured markers show the resolved observations, in which the center of mass was calculated from the 
position of Sa and Sb. The red line shows the orbit calculated by |Kohler et al. j < [2008^ . The dark blue and green lines shows two possible orbital 
fits with the extreme periods of 2.7 x 10 4 and 475 years, respectively. Right: Epochs of T Tau Sb around T Tau Sa. The coordinates are relative to 
T Tau Sa. The red line shows the orbit calculated by|Schaefer et al.|j2014|. 
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Fig. 4. The 3-cr detection limits for point-like sources near T Tau N in 
the classical imaging IRDIS mode of VLT/SPHERE and Keck/NIRC2. 
The bumps in the SPHERE data at » 250 and » 600 mas are the wings 
of the PSF from T Tau N. The Keck/NIRC2 data is obtained from the 
Keck Observatory Archive, from projects identified as N26N2 for the 
H-band and 201 OB for the K-band data. 



Fig. 5. Extended emission south of T Tau Sa. Both T Tau Sa and Sb 
were subtracted using T Tau N as a PSF reference scaled to the peak of 
T Tau Sa and Sb, respectively. The leftover artefacts are due to the Airy 
rings apparent around all three components, which cannot be totally 
eliminated because of the closeness of T Tau Sa and Sb. The leftover 
residuals have an amplitude of 2.4% and 6.6% of the peak fluxes of T 
Tau Sa and Sb, respectively. 






















